Finite thermal conductivity in ID models having zero Lyapunov exponents 
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Heat conduction in three types of ID channels are studied. The channels consist of two parallel 
walls, right triangles as scattering obstacles, and noninteracting particles. The triangles are placed 
along the walls in three different ways: (a) periodic, (b) disordered in height, and (c) disordered in 
position. The Lyapunov exponents in all three models are zero because of the flatness of triangle 
sides. It is found numerically that the temperature gradient can be formed in all three channels, 
but the Fourier heat law is observed only in two disordered ones. The results show that there might 
be no direct connection between chaos (in the sense of positive Lyapunov exponent) and the normal 
thermal conduction. 



Recent years has witnessed an increasing attention to 
the establishment of a connection between macroscopic 
phenomena such as transport coefficient and microscopic 
chaos Jl]-[|. A direct mathematical derivation has been 
proved to be very difficult, and only very simple model 
can such an approach be established Q , we have to rely 
on massive numerical simulations. There have been a 
large number of numerical works on heat conduction in 
ID systems § |TJ aim to understand what arc the nec- 
essary and sufficient conditions for a Hamiltonian system 
to obey the Fourier heat conduction law. It is found that 
an on-site potential is sufficient for a ID lattice model to 
have a finite thermal conductivity 

Albeit many progress achieved, open questions remain 
(see recent review 20 1). For example, in connecting the 
normal heat conduction with the underlying dynamics, 
some contradictions exist. On the one hand, some models 
like the ding-a-ling model [|| and the Lorentz gas model 
(with periodic and/or disordered disks) showing expo- 
nential instability, thus a positive Lyapunov exponent, 
have a normal heat conduction |Q,||,^3). On the other 
hand, the Fermi-Pasta-Ulam (FPU) model has a diver- 
gent thermal conductivity || even though it has positive 
Lyapunov exponents. Therefore, what's a role does chaos 
(in the sense of positive Lyapunov exponent) play in the 
normal heat conduction is still a unsolved problem and 
deserves further investigation. 

In this Letter, we study this problem in a series of ID 
models having zero Lyapunov exponents. Our models 
are variants of the Ehrenfest model [^J and thus called 
"Ehrenfest gas channels" . The channel consists of two 
parallel walls, a series of isosceles right triangles with 
hypotenuse along the parallel walls, and noninteracting 
particles. The two ends of the channel are put in con- 
tact with heat baths. By placing the triangles in different 
ways, we obtain different types of channels. 

The Ehrenfest model differs from the Lorentz gas 
model in underlying dynamics. The collisions of the par- 



ticles with the circles in the Lorentz gas lead to expo- 
nential separation of nearby trajectories, thus a positive 
Lyanpunov exponent, whereas collisions with the squares 
in the Ehrenfest model lead to linear separation of nearby 
trajectories, thus a zero Lyapunov exponent. 

Channel with periodic structure In this channel, the 
right triangles are placed periodically, namely, in each 
cell, we have two triangles, one on the bottom wall, the 
other on the top. The triangles are placed at the position 
of x = 1, 3, • • • (arbitrary unit). The model geometry is 
shown in Fig. 1(a). The channel of length N is iV'th 
repetition of the cell. Two heat baths with temperature 
T + and T_ are attached to the left and right ends of the 
channel, respectively . The heat bath has simple veloc- 
ity distribution Pt{ v ) = S(v — \/2T). It can be proved 
that the form of heat baths does not affect the transport 
behavior in our systems. 

To compute temperature field at a stationary state, 
we calculate time averages by dividing the configuration 
space into a set of boxes Ct |l3j . The time spent within a 
box in the jth visit is denoted by tj and the total number 
of crossings of a box Cj during the simulation is M. The 
temperature field is defined by ll^] 



T Ci = (E) Ci = 



(i) 



Then it is projected on x direction (the transport direc- 
tion). The heat flux is calculated by the change of energy 
carried through to the left and right ends by the particles, 



j = — Yae 
t M f-i ■ 



j • 



(2) 



where AEj 



(Ei, 



E n 



is the energy change at the 



jth collision with a heat bath, t jv/ is the total time spent 
for M such collisions. 
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In numerical simulation, we compute the flux for a sin- 
gle particle J\. The scaled heat flux is Jn{N) — NJ\(N) 
p3| , where N is the number of the cells. Each cell has 
length a, thus the channel has length L — Na. 
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FIG. 1. Periodic channel, (a) Geometry, (b) Temperature 
profile, (c) Heat flux J X (N) versus N. (d) ((x(t) - x(0)) 2 ) 
versus t. In (c) and (d) the bullets are the numerical data and 
the solid lines are the best fit ones. The width of the channel 
is 1.1 (arbitrary unit) and the height of right triangle is 0.6. 

In spite of the jumps at both ends the temperature 
gradient is well established and scales as dT/dx cx TV -1 
as shown in Fig. 1(b). The heat flux Ji(N) is found to 
be 



J^N) = AN C 



(3) 



with a = -1.186 ± 0.002. The thermal conductivity 
k = — J ^ul oc NJn ~ TV - 81 , which is divergent as 
one goes to the thermodynamic limit (N — > oo). 

To understand this divergent behavior, we study the 
transport property of the particles in the channel quan- 
tified by the mean square displacement ((x(t) — x(Q)) 2 ). 
An ensemble of particles (10 5 ) with the same amplitude 
of velocity (= 1) are injected into the channel in random 
directions. The best fit for the asymptotic behavior gives 
rise to 



((x(t)-x(0)) 2 )=Dt f) , 



(4) 



with = 1.672 ± 0.003 (Fig 1(d)). This means that 
the transport along x direction is neither a ballistic one 



((3 = 2) nor a diffusive one (/3 = 1). This super diffusion 
is responsible for the divergent thermal conductivity. It 
may also be the reason for the jumps near the channel 
ends in the temperature profiles. Such jumps have been 
observed in the FPU model |||?J and attributed to the 
soliton alike excitations HQ . A quantitative analysis has 
been done by Aoki and Kusnezov Jl7| more recently. 

When our model is compared with the Lorentz gas 
channel |l3| ] , it is intuitive to attribute the divergent ther- 
mal conductivity to the zero Lyapunov exponent. To 
clarify this point, we modify the channel slightly in two 
ways: (a) make the height of triangles random; (b) put 
the triangles in random position along the transport di- 
rection. The Lyapunov exponent in both variants remain 
zero because of the flatness of the triangle sides. 

Channel with right triangles of random heights The 
height of the triangle is given by 



h = ho + d * Ri, i = l,2,---,2N 



(5) 



where {Ri} are random numbers uniformly distributed 
in the interval [—1,1], d is the magnitude of disorder. 
hi < H, where H is the width of the channel. Figure 
2(a) shows the geometry of the channel. In this Letter, 
we take H = 1.1, h = 0.6, and d G [0,0.4]. 

In our calculations, the temperature and heat flux are 
averaged over 100 disorder realizations and compared 
with that one from averaged over 1000 realizations, the 
difference is found to be indistinguishable. 

Figure 2(b) shows the temperature profile for d = 0.4. 
It is a straight line with gradient dT/dx = — 0.05/7V. The 
heat flux Ji(N) is shown in Fig. 2(c). The best fit gives 
rise to a slope a = -1.992 ± 0.018. Therefore J N {N) ~ 
N^ 1 . The thermal conductivity k = -J N (N)J (dT/dx) 
is an TV independent constant, the Fourier law is thus 
justified. To see how the heat conduction changes with 
disorder, we calculate the exponent a for different values 
of d by fixing the channel length. The results are shown 
in Fig. 2(d). The bullets represent the a values obtained 
from the best fit with Eq.(|) by using N € [16,512]. 
It shows that, for a disordered channel of finite length, 
the heat conduction obeys Fourier law when the disor- 
der amplitude is large enough. In principle, in the ther- 
modynamic limit, any infinitesimal disorder will cause 
a diffusive transport, thus a normal thermal conduc- 
tion. This is demonstrated by the case with d = 0.0125, 
a = -1.724 ± 0.012 from the data N € [16,512] which 
is far from the normal thermal conduction, however, 
a = -1.999 ±0.011 from N £ [1024,32768] (the star 
in Fig. 2(d).) showing a normal thermal conductivity. 
This is similar to the mass disordered lattice model [[To] . 
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FIG. 2. Height disordered channel, (a) Channel geometry, 
(b) Temperature profile, (c) Heat flux Ji(N) versus N. (d) 
a versus d. (e) /3 versus d. d — 0.4 in (b) and (c). The solid 
lines in (d) and (e) are drawn to guide eyes. 

We compute ((x(t) — x(0)) 2 ) and find that for all values 
of disorder, it can be best fitted by Dt" asymptotically. 
f3 as a function of the disorder d is plotted in Fig. 2(e). 
It is seen that for any finite value of d, the slope (3 is very 
close to unity, which means that the particles moves dif- 
fusively in the channel, thus the heat conduction in this 
channel obeys Fourier law. 

Thermal conductivity k versus temperature To = 
(T+ + T_)/2 is plotted in Fig. 3(a). It is found that 
k ~ T 7 , and the best fit gives rise to 7 = 0.501 ± 0.002. 
The normalized temperature profile T*(x) — T(x)/Tq 
is shown in Fig. 3(b), which indicate that dT/dx = 
-0.02T /L. 

Channel with triangles at random positions The posi- 
tion of the triangle is made random, namely, Xi = d* Ri, 
where x^ is the position away from the periodic structure 
shown in Fig. 1(a). Figure 4(a) is the schematic illus- 
tration of the geometry. The numerical simulations show 
that the temperature gradient is well established and is 
similar to Fig. 2(b). The heat flux J\{N) is described 
by Eq. (|[). We plot the exponent a versus d in Fig. 
4(b). It tells us the trend to normal thermal conductiv- 
ity (a = —2) in a long length limit (star) and a large 
disorder limit (bullets). 

As an independent check, we also calculate the integral 
of the current-current correlation function in the Green- 



Kubo formula. The integral is found to be convergent in 
cases with disorder but divergent in the case with peri- 
odic geometry shown in Fig. 1 (a) . 




FIG. 3. Height disordered channel, (a) Thermal conduc- 
tivity k versus temperature T = (T+ + T_)/2. (b) The 
normalized temperature profile (T* = T /To) for six differ- 
ent temperature scales. To = 0.01,0.1,1,10,100, and 1000, 
respectively. Disorder d = 0.4. 

Theoretical analysis Suppose the path length distri- 
bution of particles from left to right (or vice versa) is 
/l{1) in a channel of length L, namely, there are Sn par- 
ticles whose path length lies in the interval [1,1 + dl\. 
bn/n — /i(Z)dL is determined merely by the struc- 

ture and the length of the channel. Two heat baths of 
temperature Tj, and Tr are put to the left and right 
ends, respectively. In a time period of t there are n 
particles exchanged between two heat baths. The to- 
tal time spent to reach the right heat bath from the left 
one is t LR = n(l) J °° ±P Tl (v)dv, where (I) = J °° lf L (l)dl 
is the average path length from the left heat bath to 
the right one, and Pt{v) the velocity distribution func- 
tion of heat bath at temperature T . Similarly, the to- 
tal time of n particles from the right bath to the left 
one is tut = n(l) J °° ^Pt 2 {v)dv. The total time is 
t = tLR + iijL- The energy exchange between two heat 
baths is E = n J °° ^{P Tl {v) - P T2 {v))dv = n{T x - T 2 ), 
the heat flux for the channel of length L per particle is 
thus given by: 



J l( L) = f 



Ti -T 2 



(6) 



(I) f o °° l(P Tl (v) + P T2 (v))dv- 

From Eq(^) , it can be seen that whether the heat con- 
duction obeys the Fourier law or not does not depend on 
the types of heat bath, it depends only on (/) - the trans- 
port property. For instance, if the system is diffusive, 
then (1) oc L 2 and the heat flux Ji(L) cx L~ 2 . This is 
what we see in numerical calculations (Fig. 2(c)). More- 
over, for a given geometry, i.e. (I) is determined, the heat 
flux and heat conductivity are determined by the prop- 
erty of the heat baths. If we change the temperature 
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of heat baths q times, then for the simple heat bath we 
used and the Gaussian heat bath, it can be shown that, 
the heat flux changes q 3 / 2 times. Because the tempera- 
ture gradient dT/dx = Const. To (see Fig. 3(b)), thus the 

1/2 

thermal conductivity k changes with temperature T , 
this agrees with our numerical finding in Fig. 3(a). 
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FIG. 4. Position disordered channel, (a) Geometry, (b) a 
versus d. {xi} — d * Ri, d is disorder. The solid line in (b) 
is drawn to guide eyes. In (b) the bullets are the data from 
best fit by using N G [32, 1024], while "*" at d = 0.1 is the 
one by using N € [8192, 32768]. 

In summary, we have studied heat conduction in three 
different ID Ehrcnfcst channels. The temperature gradi- 
ent can be formed in all cases. However, a finite thermal 
conductivity can be reached only when the disorder (ei- 
ther in position or in height) exists. As the Lyapunov 
exponents are zero in our model, we thus conclude that 
the finite thermal conductivity might have nothing to do 
with the underlying dynamics. Most recent study on heat 
conduction in channels with irrational triangles supports 
this argument p2[ . 
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